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ABSTRACT 



Aims. We intend to assess the most plausible scenarios for generating large amounts of dust in high-z quasars (QSOs) on the basis of 
observationally derived physical properties of QSOs at z > 6. 

Methods. We use a chemical evolution model to compute the temporal progression of quantities such as the amount of dust and gas, 
stellar masses, star formation rates (SFRs) and the metallicity for various combinations of the initial mass function (IMF), the mass 
of the galaxy, dust production efficiencies, and the degree of dust destruction in the ISM. We investigate the influence of the SFR on 
the evolution of these quantities, and determine the earliest epochs at which agreement with observations can be achieved. We apply 
the obtained results to individual QSOs at z > 6. 

Results. We find that large quantities of dust can be generated rapidly as early as 30 Myr after the onset of the starburst when the 
SFR of the starburst is > 10 3 M Q yr -1 . The amount of dust and several other physical quantities of individual QSOs at z > 6 are 
satisfactorily reproduced by models at epochs 30, 70, 100, and 170 Myr for galaxies with initial gas masses of 1-3 X 10" M Q . The 
best agreement with observations is obtained with top-heavy IMFs. A sizable dust contribution from supernovae (SNe) is however 
required, while at these epochs dust production by asymptotic giant branch (AGB) stars is negligible. Moderate dust destruction in 
the ISM can be accommodated. 

Key words, galaxies: high-redshift - galaxies: starburst - galaxies: evolution - ISM: evolution - quasars: general - stars: massive 



1. Introduction 

Studying QSOs and their host galaxies at high redshift (z > 6) 
is important to gain deeper insight into the formation and evolu- 
tion of galaxies, the origin of dust production, and the build up 
of stellar bulge masses in coevolution with supermassive black 
holes (SMBHs). 

While the m ost distant known QSO, Jl 14816.64+525150.3 
(iFan et al.l 12001 herafter Jl 148+5251), is at z = 6. 4, several 
tens of QSOs have been discovered at z ~ 6 (e .g., IFan et al.l 
12004 120061: IWillott et al.ll2007L Uiang et al.ll20Toh . Most of the 
observed QSOs at this redshift, where the epoch of cosmic 
evolution is ~ 1 Gyr, exhibit extreme physical properties such 
as very high far-i nfrared (FIR) luminosities which imply large 
dust masses (e.g.. lOmont et al.ll200lL [20031: ICarilli et alJl200lt 



distributed over kilo -parsec scales throughout the host galaxy 
dLeipski et al . 2010). The amount of cold dust inferred is about 
a few times 10 8 M (e.g.. iBertoldi et al.ll2003at iRobson et al.l 



I2004t iBeelen et al.ll2006f iMichalowski et all 1201 0l) . The domi 
nant source of the high FIR luminosity is believed to be dust 
heate d by intense star formation in the circumnuclear regio n 
(e.g. JCarim et al.ll2004tlRiechers et al.ll2007tlWang et"aT1l2 008). 
Detec tion of [Cll] line emission at 158 /mi (Maiolino et al. 
12001 within a central region with radius ~ 750 pc of the host 
galaxy of Jl 148+5251 also implies a high star formation rat e 
surface density of 1000 M yr 1 kpc 2 dWalter et all [2009). 
IWang et al.l $201$) derived SFRs between 530-2300 M Q yr" 1 
from observations of a sample of QSOs at redshift z > 5. 
Observations of s t rong metal emission of high-z QSO s (e.g. , 
Bart h et all 120031. iDietrich et al.l l2003t IMaiolino et all 120031; 
iBecker et al.l 12006) indicate strong star forming activ ity in the 
QSO hosts and solar or sup e rsolar metall i city (e .g., IFan et al.l 
2003; iFreudling et ail 120031; Uuarez et all 120091) . Theoretical 
studies of the gas metallicity of QSO h osts also predict supers o- 
lar metallicities for z = 5-6 QSOs (e.g. JPi Matteo et al.l l2004). 



I Bertoldi & Coxll2002l). and SMBHs with masses > 10 9 M H (e.g., 
Barth et al.ll2003t IWillott et al.ll2003HVestergaardll2004 . 

Observations of QSOs have shown that dust emission at 
near-infrared (NIR) wavelengths arise from warm a nd hot dust 
(T <, 1000 K) assemble d within a few parsec (e.g. jHines et al.l 
2006; Uiang et all 120061) . The NIR emission is believed to be 
powered by th e active galactic nu cleus and related to the QSO 
activity (e.g.. iPolletta et al.l |2006T) . However, two QSOs at z ~ 
6 without detectable em ission from hot dust have been found 
dJiang et al.ll2006l l2010h . It has been proposed that these QSOs 
are at a too early evolutionary stage to have built up signifi- 
cant amounts of hot dust. Alternative scenarios including for 
example the destruction of the hot dust or dust misalignmen ts 
from the SMBH h ave also been discussed (lHao et al.l l2010allbl : 
iGuedes et aTll2010h . 

The FIR luminosity of Lhr ~ 10 12 ~ 13 L H is attributed to cold 
dust (T ~ 30-60 K) (e.g., Wang et al. 2008) which is probably 



The high inferred SFRs imply short timescales (< 10 yr) 
of the starburst (e. g. IBertoldi et al.l l2003aL IWalter et all [2004; 



iDwek et all 120071; iRiechers et al.1 120091) . and consequently a 
young age of the QSOs. An early evolutionary stage of z > 4 
QSOs has also been suggested from studi es of extinction curves 
of broad absorption line QSOs (e.g., iGallerani et al.l 1201 Ol) 
which turned out t o be best fitted with extinction curves for 
SN-like dust (e.g.. IMaiolino etal.ll2004l 120061: IGallerani et all 
I2010I) . This suggests SNe as the preferential source of dust 
at early epoch s (e.g. , IDwek 1998t Morgan & Edmundsl 12 003; 
iHirashita et all 120051: IDwek et alJ I2007L IDwek & Cherchnefl 
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2010), even though the dust productivity of SNe is poorly con- 
strained (for a review see Gall et al. in prep). Th e dust in high- 
z QSOs could also be grown in the ISM (e.g., |Drainej[2009; 
iMichalowski et all l2Qlof iPipino et aD 1201 ll) . Finally, a domi- 
nant dust production by asym ptotic giant branch stars has been 
claimed jValiante et alj2009t) . 

Molecular gas masses of the order of ~ 1-2.5 x 10 10 M 
have been inferred from detections of high excitation CO line 
emiss ion in QSOs at z > 5 within a ~ 2.5 kpc radius region 
(e ^lBertoldieTaIll2003bt IWalter et alJl200l 12001: IWangetalJ 
2010). The dynamical masses inferred from these CO observa- 
tions are a few times ~ 10 10-11 M which sets an upper limit on 
stellar bulge masses. These however are roughly two orders of 
magnitude lower th an required from the pr esent day black hole- 
bulge relation (e.g jMarconi & Huntl2003l) . It therefore has been 
proposed that the formation of the SMBH occurs prior to the 
formation of the stellar bulge. QSOs will then have to accrete 
additional materia l to build up the required bulge mass (e.g., 
IWalter et alJl200l: iRiechers et al1l2009t IWang et al.ll20Tol) . For 
QSOs at z > 6 super-Eddington growth on timescales shorter 
than ~10 8 yr seem to be requi red to form a SMBH > 10 9 M 
(e.g., iKawakatu & W ada 2009]). It has also been predicted that 
QSOs at z ~ 6 likely have formed in dark matter h alos of 10 12 ~ 13 
M Q (e.g. JLi et alfeOOllKawakatu & Wadall2009l) . 

In lGall et alJ d201 ll herafter fpaper J) we developed a chemi- 
cal evolution model to elucidate the conditions required for gen- 
erating large dust masses in high-z starburst galaxies. We showed 
that galaxies with masses of 1-5 x 10" M are suitable for en- 
abling the production of large amounts of dust within ~ 400 Myr. 
In the present paper we apply this model to QSOs at z > 6. We 
perform more detailed comparison between model results and 
values inferred from observations of z > 6 QSOs to identify the 
most likely scenario. Furthermore, we consult additional param- 
eters such as the H2 mass and the CO conversion factor for more 
refined evaluati ons. In particular, calculations with higher SFRs 
than in lPaper j are considered. We aim to determine the earliest 
epochs at which the model results are in agreement with those 
from observations. 

The structure of the paper is as fol lows: In Sec.|2]we briefly 
review the model developed in Paper I. A detailed analysis of the 
results is presented in Sec.[3]followed by a discussion in Sec. [4] 



2. The model 

The galactic chemical evolution model from Paper I is self- 
consistent, numerically solved and has been developed to as- 
certain the temporal progression of dust, gas, metals, and di- 
verse physical properties of starburst galaxies. The incorporated 
stellar sources are AGB stars in the mass range 3-8 M and 
SNe. A differentiation between diverse SN subtypes has been 
implemented. Their roles as sources of dust production, dust de- 
struction or suppliers of gas and heavy elements are taken into 
account. The lifetime dependent yield injection by the stellar 
sources, as well as dust destruction in the ISM due to SN shocks 
are also taken into account. Moreover, the formation of a SMBH 
is considered. Due to the very high SFRs of the starbursts, infall 
of neutral gas will only effect the system for comparable high 
infall rates. Thus, gas infall and outflows are neglect ed. Possible 
caveats of such an approach are discussed in Paper I. The model 
allows investigations of a broad range of physical properties of 
galaxies. 

The prime parameters are summarized in the following. 



Table 1. Model parameters 



Parameters 


Value 


Unit 


Description 


<Aini 


3 x 10 3 , 1 x 10 4 


M G yr" 1 


Star formation rate 


^SMBH 


3 x 10" 


M Q 


Mass of the SMBH 


fSMBH 


1 x 10 s 


yr 


Growth timscale 








for the SMBH 



- Three different possible prescriptions for the stellar yields 
of SNe are implemente d, i.e., (i) stellar evolution models 
bv lEldridge etal] d2008| ) (referred to as 'EIT 08M'), (ii) ro- 
tating stellar models bv lGeorgy et al.l d2009l). or (i ii) nu cle- 
osynthesis models b y either IWooslev & Weaver! d 19951) or 
Nom otoet al.1 d2006|). The stellar yields for AG B stars are 
taken from lvan den Hoek & Groenewegenldl997l) . 

- We differentiat e between five different IMFs. These are a 
iSalpeteri (fl95l IMF, a top-he avy, and a mass-heavy IMF, 
as well as IMFs dLarsoJ l998) with characteristic masses of 
either m c h = 0.35 (Larson 1) or ;« c h = 10 (Larson 2). 

- The SFR at a ce rtain epoch is given by the Kennicutt law 
dKennicuttll 19981) as ifr(f) = tA ini (M ISM (t)/M ini ) k , where ifr^ 
is the initial SFR of the starburst, Mism(0 is the initial gas 
mass of the galaxy and k = 1.5. 

- The amount of dust produced by SNe and AGB stars is cal- 
culated using the dust formation efficiencies discussed in 
iPaper 1 For SNe three different dust production efficiency 
limits are determined, i.e. a 'maximum' SN efficiency, a 
'high' SN efficiency, and a 'low' SN efficiency. The 'maxi- 
mum' SN efficiency originates from theoretical SN dust for- 
mation models, and corresponds to dust masses of approx- 
imately 3-10 x 10 -1 M . Similar dust masse s have been 
obser ved in SN remn ants such as Cas A (e.g., iDunne et al.l 
2009) or Kepler (e.g.. lGomez et al.ll2009t) . Dust destruction 
in reverse shock interaction of about 93 % has been applied 
to the 'maximum' SN efficiency, to obtain the 'high' SN ef- 
ficiency. The amount of dust for instance is ~ 2-6 x 10~ 2 
M , which is also com parable to some observations of older 
SN remnants dPaper ll see references therein). The 'low' SN 
efficiency is based on SN dust yields (on average about 3 x 
10~ 3 M ) inferred from observations of SN ejecta. 

- Dust destruction in the ISM is implemented in terms of the 
mass of ISM material, M \, swept up by a single SN shock 
and cleared of the containing dust. 

For calculations in this pa per most parameters have the same 
settings as defined in Paper I. We apply the models where the 
formation of a SMBH has been included. A constant growth rate 
has been estimated based on the final mass of the SMBH and the 
considered growth timescale. In this paper the SMBH growth is 
considered with a shorter growth timescale and calculations are 
performed with higher initial SFRs. For the SN yields we only 
consider the case of EIT08M. The parameters which differ from 
those used in lPaper B are listed in Table [T] 

3. Results 

In this section we present the results of models calculated within 
short timescales after the starburst. 

A short enrichment timescale of a few times 10 7 yr for an 
intense starburst with a SFR of ~ 3 X 1 3 M yr" 1 has bee n 
proposed by e.g iBertoldi et all d2003al) . IWalter et al] (|2004), 
iDwek et alJd2007l)jRiechers et al.ld2009l) Owing to this sugges- 
tion we are interested in whether the observed large dust masses 
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Fig. 1. Relation between dust mass and stellar mass at an epoch of 30 Myr, for various initial gas masses and IMFs. Calculations 
are performed for a 'maximum' SN efficiency and dust destruction in the ISM with M c i = 100 M (left panel) and M c \ = (right 
panel). The colored symbols are obtained for different initial gas masses, M m \, SFRs, and IMFs. The size of the symbols is scaled 
by M ini . Calculations are made for M ini = 1.3 x 10 12 M G (largest symbol), M ini = 5 x 10 11 M , M ini = 3 x 10 11 M G , M ini = 1 x 10 11 
M , and M- m \ = 5 x 10 10 M (smallest symbol). The crosses correspond to calculations for a initial SFR ij/ m i = 10 3 M yr _1 , the filled 
circles to t^ini = 3 x 10 3 M yr~\ and the stars to ij/ m \ = 10 4 M yr _1 . The black, green, cyan, magenta, and blue colors denote the 
Salpeter, mass-heavy, top-heavy, Larson 1, and Larson 2 IMF, respectively. The dark grey region indicates the mass range of stellar 
masses and dust masses derived from observations of QSOs at z > 6. The vertical dashed lines represent the lower and upper limits 
of the observed stellar masses. The light grey area illustrates the whole mass ranges derived from observations of QSOs > 5 and 
accounts for uncertainties in the derived quantities. The horizontal dashed lines mark the lower and upper mass limit of the derived 
dust masses. 



in excess of 10 M can be reached within about 100 Myr. 
Consequently we performed calculations with an initial SFR for 
the starburst with tfr m [ = 3 x 10 3 M yr~' for galaxies with initial 
gas masses M ini = 5 x 10 10 M , M ini = 1 x 10 11 M Q , M ini = 3 x 
10 11 M Q , and M ini = 5 x 10 11 M Q . For the most massive system 
with M ini = 1.3 x 10 12 M an initial SFR ifr ini = 1 x 10 4 M yr _1 
is adopted. We included the results for a low er initial SFR of 10 3 
M yr~' from models computed in Paperl for comparison. 

In Pa per il we analyzed the evolution of the amount of dust 
and various physical properties, and found that these are strongly 
dependent on the mass of the galaxy. Moreover, for a given ini- 
tial SFR all quantities evolve faster in less massive galaxies. In 
this paper we perform detailed comparisons between calculated 
and observed values of the total dust mass, Mj, the stellar mass, 
M*, the SFR, iff, and the metallicity, Z. We identified the shortest 
epoch, where some model results are in accordance with obser- 
vations to be 30 Myr. Furthermore, we discuss quantities such 
as the CO conversion factor, the gas-to-H2 mass ratio, and the 
possible amount of molecular hydrogen. 

3. 1. Dust and stellar mass 

In Fig. Q] we present the results for the mass of dust versus the 
stellar mass for galaxies with different initial gas masses and 
initial SFRs at an epoch of 30 Myr. The displayed models are 
computed for a 'maximum' SN efficiency. Dust destruction in 
the ISM is considered for values of M c] = 100 M (left panel) 
and M c i = (right panel). 

The dark grey region represents the mass ranges of the stellar 
mass and dust mass derived from observations of QSOs at z > 6. 
The lower and upper limits of the stellar mass are estimated by 
subtracting the molecular gas masses, Mh, from the total dynam- 
ical masses, M^n- Valu es for Md yn and Mh, are based on data 
from [Wang et all d2010l and references therein) for three QSOs 
at z > 6. For an estimation of Md yn an inclina tion angle i = 65° o f 
the gas disk is taken for QSO Jl 148+5251 (IWalter et al.ll2004h . 
while i = 40° similar to IWang et alj d2010h is applied to the 
remaining two QSOs. We adopt the lower and upper limits for 



the du st masses from Beel en et al ] (l2006h and Michalowski et al. 
d2010 ). The light grey region covers the range of derived stel- 
lar masses and dust masses from observa tions of QSOs > 5 
dWang et al.ll2010l:lMichalowski et all2 0101. The boundaries for 
the stellar masses are estimated similar to the QSOs at z > 6 
(with i = 40° for deriving Md yn ). We set the lower dust limit to 
10 8 M to account for the uncertainties of derived dust masses 
from observations. 

Despite the short time span of 30 Myr, it is evident that most 
models are within the plausible mass ranges illustrated by the 
light and dark grey regions. This signifies a rapid build-up of a 
large amount of dust, provided SNe produced dust with a 'maxi- 
mum' SN efficiency. For galaxies with M; n i = 1-5 x 10 11 M all 
models with an initial SFR of 3 x 10 3 M yr are in agreement 
with the observed values for the stellar masses for QSOs at z > 
6. The requirements for Mj are best accomplished with either a 
top-heavy, mass-heavy or Larson 1 IMF for both values of M c \, 
In a galaxy with M;„i = 1 x 10 11 M Q the amount of dust reached 
with a Larson 2 IMF and M c i =100 M also matches with the 
dark grey region. Models for either a 'high' or 'low' SN effi- 
ciency did not reach 10 8 M Q of dust. Only in the most massive 
galaxy (Mi„i = 1.3 x 10 12 M ) and for top-heavy IMFs with a 
'high' SN efficiency an amount of dust > 10 8 M Q is obtained. 

In Fig.|2]we illustrate the results for dust and stellar masses at 
an epoch of 100 Myr. We present models for a 'maximum' SN 
efficiency (top row) and a 'high' SN efficiency (bottom row), 
while dust destruction in the ISM is considered for a M c i = 800 
M Q (left column), M c i = 100 M (middle column), and M d = 
(right column). We carried out calculations for a 'low' SN effi- 
ciency, but the obtained dust masses of these models remained 
below 10 s M . 

At these early epochs the stellar mass, M*, is higher for mod- 
els with an initially larger SFR (at fixed IMF and M;„i). The stel- 
lar mass is also larger for IMFs biased towards low mass stars 
(at fixed M;,,; and i/'ini). It is interesting to note that in the less 
massive galaxies (0.5-1 x 10 M ) dust masses obtained for 
the higher initial SFR (i// iai = 3 x 10 3 M yr~') are lower than 
dust masses obtained for the lower SFR (i/fj m = 10 3 M Q yr' 1 ). 
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Fig. 2. Relation between dust mass and stellar mass at an epoch of 100 Myr for various initial gas masses and IMFs. Calculations 
are performed for a 'maximum' SN efficiency (top row) and a 'high' SN efficiency (bottom row). Dust destruction in the ISM is 
considered for a M c i = 800 M (left column), M c \ = 100 M (middle column), and M c \ = (right column). The colored symbols are 
obtained for different initial gas masses, Mj n ;, SFRs, and IMFs. The size of the symbols is scaled by Mj n ;. Calculations are made for 
Mi n i = 1.3 x 10 12 M Q (largest symbol), M ini = 5x10" M , M ini = 3 x 10 11 M Q , M ini = 1 x 10 11 M G and M ini = 5 x 10 10 M (smallest 
symbol). The crosses correspond to calculations for a initial SFR \fi^ = 10 3 M yr~', the filled circles to \j/ m i = 3 x 10 3 M yr~', 
and the stars to ifr m i = 10 4 M yr _1 . The black, green, cyan, magenta, and blue colors denote the Salpeter, mass-heavy, top-heavy, 
Larson 1, and Larson 2 IMF, respectively. The dark grey region indicates the mass range of stellar masses and dust masses derived 
from observations of QSOs at z > 6. The vertical dashed lines represent the lower and upper limits of the observed stellar masses. 
The light grey area illustrates the whole mass ranges derived from observations of QSOs > 5 and accounts for uncertainties in the 
derived quantities. The horizontal dashed lines mark the lower and upper mass limit of the derived dust masses. 



Moreover, in these galaxies the amount of dust reached at an 
epoch of 30 Myr (see Fig.[fl and for M cl = 100-800 M is also 
higher than that seen at the epoch of 100 Myr for same M c \. 

We find that the stellar masses for models with an initial SFR 
•Ami = 1 _ 3 x 10 3 M yr _1 are within the observed region for 
z > 5 QSOs. For some models with iff ini = 3 x 10 3 M yr -1 , 
stellar masses are within the mass range for z > 6 QSOs. This in 
particular applies to systems with either M m \ = 0.5-1 x 10 1 1 M 
(all IMFs) or M ini = 3-5 x 10 11 M with top heavy IMFs. Stellar 
masses within the dark grey area are also found with i/^i = 10 3 
M yr~' for galaxies with either M mi = 3-13 x 10 11 M and top 
heavy IMFs or for the less massive galaxies in combination with 
IMFs favoring low mass stars. 

In the case of M c \ = 800 M and for a 'maximum' SN effi- 
ciency most models with M mi = 3-13 x 10 11 M and if/ m \ = 10 3 
M yr~' fit within the dark grey region. However for the higher 
initial SFR Mj is within or close to this zone only for galaxies 
with M ini = 3-5 x 10 11 M and top-heavy IMFs. For M d = 100 
M and a 'maximum' SN efficiency the dust mass obtained in 
a galaxy with M in i = 1 x 10 11 M , tff ini = 3 x 10 3 M yr _I and 
for top-heavy IMFs is in agreement with observations, while the 
dust masses in the more massive galaxies for some IMFs and 
SFRs are higher than required. In the case of no dust destruction 
the dust masses reached for some IMFs and SFRs are able to 
match within the dark grey area also in the least massive galaxy. 

We find that in case of a 'high' SN efficiency and for i(/ m i = 
3 x 10 3 M yr~' in galaxies with initial masses 3-5 x 10 11 M 



and top-heavy IMFs high dust masses are possible, even if dust 
destruction is included (i.e., M c i = 0-100 M ). 

3.2. Metallicity and SFR 

We next present the obtained metallicities and SFRs at the time 
of observation for the models discussed above. 

Fig. [3] depicts the metallicity versus SFR at epochs of 30 
Myr (left panel) and 100 Myr (right panel). With respect to 
observations of QSOs > (5) 6 we marked the range of de- 
rived values as a dark grey shaded zone. The lo wer and upper 
limits o f the SFR are ba s ed on observations by Bertoldi et al.l 
d2003al) and IWang et all d2010l) . We set the lower limit for 
the metallicity at the solar value and the upper limit at 5 Z . 
This is based on the i nferred solar or supersolar metallicities 
in high-z QSOs (e.g.. iBarth et alj 120031; IPietrich et all [2003; 
| Fan et alJ 12003) iFreudling et alJ l2003t iMaiolino et all 12003; 
iDi Matteo et alJ l2004t iBecker et alJ 120061: lJuarez et alJ 12009). 
We note that there are no strong constraints on the upper limit 
and therefore the zone above 5 Z is marked as light grey shaded 
region to account for the uncertainty. 

We find that at an epoch of 30 Myr high metallicities in the 
less massive galaxies are already reached. The best result is at- 
tained by a system with M;„i = 1 x 10 11 M , ip mi = 3 x 10 3 M 
yr -1 , and IMFs biased towards higher masses. For a galaxy with 
Mi n i = 5 x 10 10 M all models with either the same ifr m i or with 
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Fig. 3. Relation between metallicity and SFR at epochs of 30 
Myr (left panel) and 100 Myr (right panel). The colored sym- 
bols are obtained for different initial gas masses, Mi,,;, SFRs, and 
IMFs. The size of the symbols is scaled by M m \. Calculations are 
made for M ini = 1.3 x 10 12 M Q (largest symbol), M ini = 5 x 10 11 
M Q , MM = 3 x 10 11 M , Mi n i = 1 x 10 11 Mq, and M ini = 5 x 
10 10 M Q (smallest symbol). The crosses correspond to calcula- 
tions for a initial SFR i^; n i = 10 3 M yr , the filled circles to i^i n ; 
= 3 x 10 3 M yr" 1 and the stars to ifr^ = 10 4 M yr 1 . The black, 
green, cyan, magenta, and blue colors denote the Salpeter, mass- 
heavy, top-heavy, Larson 1, and Larson 2 IMF, respectively. The 
dark grey shaded region indicates the range of the metallicity 
and SFR based on observations of QSOs at z > 6. The vertical 
dashed lines represent the lower and upper limits of the observa- 
tionally derived SFRs. The light grey shaded area accounts for 
the uncertainty of the upper limit of the metallicity. The horizon- 
tal dashed lines mark the lower and possibly upper limit of the 
metallicity. 



the lower initial SFR, and top-heavy IMFs are within the dark 
grey shaded region as well. 

At an epoch of 100 Myr the metallicity has increased in all 
models, while the SFR in the less massive galaxies has signif- 
icantly decreased. The models for M;„j = 3-5 x 10 11 M Q , 0mi 
= 3 x 10 3 M yr -1 , and top heavy IMFs constitute the best re- 
sults. In galaxies with M in ; = 3 x 10 11 M Q , the same initial SFR, 
and either a mass-heavy or Larson 1 IMF the obtained values for 
Z and ifr(t) are also in agreement with the observed values. The 
metallicities in the low mass galaxies which give the best agree- 
ment at 30 Myr are now shifted above the upper limit, while the 
SFRs remain in the observed range. The models for a galaxy 
with Mini = 1 x 10 11 M , a lower initial SFR of 10 3 M yr 1 , 
and top-heavy IMFs at this epoch (100 Myr) reach sufficiently 
high metallicities, while high enough SFRs are sustained. 

3.3. CO conversion factor and gas-to-H 2 mass ratio 

To evaluate the calculated models, we additionally consider the 
relation between the gas-to-H2 mass ratio and the CO conversion 
factor used to derive the molecular gas mass in a galaxy. 

Detections of high excitation CO line emission in QSOs at z 
> (5) 6 indicate the presence of 0.7-2 . 5 X 10 10 M of mo l ecular 
hydrogen (e.g.. IBertoldi et alj|2003bt IWalter et alj |2003l 12004 
iRiechers et al.l2009tlWang et al.l2 010). This molecular gas mass 




Fig. 4. CO conversion factor versus gas-to-IL ratio at epochs 30 
Myr (top panel) and 100 Myr (bottom panel). The solid lines cor- 
respond to calculations of a as a function of the gas-to-FL ratio 
?7g,H2 for a CO line luminosity of ^co(i-O) = ^-7 X 10 10 K km s _1 
pc 2 . Calculations are performed for different IMFs and galax- 
ies for a range of different initial gas masses Mi,,;. The thickness 
of the lines is scaled by M; n i as indicated in the upper panel. 
The black and cyan colors denote the Salpeter and top-heavy 
IMF, respectively. The arrow indicates the shift of a for calcu- 
lations with the lower L' C0 ^ l Q) = 1.5 x 10 ll) K km s _I pc 2 , and 
A a is the difference of a between the higher and lower L' co ^ l 0) . 

Calculations are shown for models with i^; n i = 3 x 10 3 M yr 1 , 
except for the model for the most massive galaxy for which ifr^ 
= 10 4 M yr . The grey shaded region signifies the possible 
range of a and 77 g ,H 2 ■ The horizontal black dashed lines mark the 
values of a = 0.8, 1 and 4.6 M (K km s _I pc 2 ) -1 . 



is derived from the relation Mh, = a x £^0(1-0) ' wnere a i s tne 
conversion factor between the low excitation CO J = 1-0 line 



luminosity L' c 



'CO(l-O) 



and Mh 2 ■ For spiral galaxies a is typically 



4.6 M (Kkms- 1 pc 2 ) 1 (e.g.. lSolomon & BarrettlfT99l . while 
for the centre of nearby ultra luminous starburst galaxies a con- 
version f actor of a = 0.8-1 M (K km s _1 pc 2 ) 1 is appropri- 
ate (e.g., iDownes &Solomod fl998). The latter value of a is 
usually used for e. g., high-z QSOs ( e.g., IBertoldi et al.l l2003b; 
IWalter et alJ 12003: Wan g etalJl201 Oh Ultra Luminous Infrared 
Galaxies (ULIR Gs) dYan et al] 120101) or for high- z sub-mm 
galaxies (SMGs) dTecza et al.l2004l:lGreve et al.l200 5). However 
a is not well known in the case of very high excitation. 
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In our models we have computed the total (H + He) gas 
mass Mq which remains in the galaxies at a given epoch. The 
molecular gas mass, Mh 2 , constitutes a certain fraction of the to- 
tal gas mass, Mq. Hence we introduce the gas-to-H2 mass ratio 
as 77g,Hi = Mg/Mu 2 . The CO conversion factor can thereby be 
expressed as a function of rj g ,h, as 



a = 



Mq 



(1) 



'CO(l-O) 



where 77 gi H 2 > 1 is kept as a free parameter. In ULIRGs and 
SMGs a major fraction of th e gas is believed to exi st in form 
of molecular hydrogen (e.g. JSanders & Mirabellfl996l) . For ex- 
ample a value for of ~ 1 has been found f or the z = 3 radio 
galaxy B3 J2330+3927 dDe Breuck et al.ll2003l) . This might also 
be the case for QSOs and suggests a gas-to-H2 ratio between 1 
and 2. 

In Fig. |4] we show the results for a as a function of r]^ 2 
with tfr mi = 3 x 10 3 M yr~' for models with M ini < 5 x 10 1 ' 
M and with i/fj n ; = 10 4 M G yr _1 for the most massive galaxy. 
Calculations are performed for two different epochs; 30 Myr (top 
panel) and 100 Myr (bottom panel). The IMFs involved are the 
top-heavy IMF and the Salpeter IMF. We adopt a CO line lu- 
minosity L' CO(l 0) = 2.7 x 10 10 K km s _1 pc 2 which is based 
on the derived values of Jl 148+5251 and J0840+5624 (e.g., 
iBertoldi et al]l2003bl: I Walter et aUl2003l:IWang et alJl201flh . 



The difference of a from calculations with a lower L' n , 
= 1.5 x 10 10 K km s- 1 — 2 ^ 



^CO(l-O) 

(i.e., L' CO(l l)) = 1.5 x 10 1U K km s _1 pc 2 ) is indicated by the 
arrow in Fig. [4] The grey shaded area signifies a possible range 
for a and Tj gz a 2 as discussed above. 

For a fixed value of a the gas-to-H2 ratio increases with in- 
creasing initial mass of the galaxy. This is as a consequence of 
the larger amounts of gas mass remaining in the more massive 
galaxies at the epochs of interest (see also Paper I). Conversely, 
for a fixed ?7 g ,H 2 , a increases with increasing M m \. The maximum 
value of a is obtained for = 1, i.e., Mq = Mh, . We find that 
at both epochs, the maximum value of a for the less massive 
galaxies is lower than ~ 4.6 M (K km s pc 2 ) -1 . For a given 
Mi n i, a, and T] g ,H 2 are lower at later epochs. For a lower £ CO(1 _ 0) , 
a shifts to higher values for a given 7/ g H 2 ■ 

At an epoch of 30 Myr the values for a and r] gt n 2 are similar 
for all IMFs and galaxies with M ini > 1 x 10 11 M Q , while the 
difference becomes larger with decreasing M m \. Feasible values 
of a and ?7g,H, are possible for galaxies with Mj n ; = 1 x 10" 
M and the higher value of L' CO(l _ 0) . For top-heavy IMFs rj gz a 2 

= 1 results in a maximum a of ~ 2.3 M (K km s pc 2 ) -1 , 
while for a = 0.8 M (K km s pc 2 ) -1 , the fraction of molecular 
hydrogen is about one third of the total gas mass. In the least 
massive galaxy (Mj n ; = 5 x 10 ll) M ) and for a top-heavy IMF a 
~ 0.8 M (K km s pc 2 ) -1 presupposes that all the gas in this 
system is in the form of molecular hydrogen. In more massive 
systems with M ilu = 1-3x10" M , a value of a x 0.8-1 M (K 
km s _1 pc 2 ) -1 presumes that the molecular hydrogen constitutes 
only a small fraction of about 1/10-1/20 of the total gas mass. 

At an epoch of 100 Myr a clear separation between the 
IMFs is noticeable. For a Salpeter IMF the galaxies underwent a 
stronger gas exhaustion than for a top-heavy IMF, which is more 
significant for the less massive galaxies. As for the epoch at 30 
Myr the system with M ini = 1 x 10 11 M and top-heavy IMF 
is plausible , i.e., for a ~ 0.8 M (K km s _1 pc 2 ) -1 the gas-to- 
H 2 ratio T] g ^ 2 = 2. For the galaxies with M mi = 3-5 x 10 11 M 
and top-heavy IMF we obtain a = 1.4-1.5 for a corresponding 
gas-to-H2 ratio rj g ii 2 = 5-10, resulting in a molecular mass of 



Mh 2 ~ 3.7 x 10 M . Alternatively, a higher value for a up to 
4.6 results in a lower T] gt n 2 = 2-4. It is noteworthy that for the as- 
sumed Lco(i-O) = 2 - 7 x 10 10 Kkms _1 pc 2 , a = 4.6 M (Kkm s -1 
pc 2 ) 1 implies M H2 = 1.2 x 10 11 M . The likelihood that such a 
high Mh, could have been built up within a short timescale of 
30-100 Myr however is unclear. 

4. Discussion 

4. 1. Individual QSOs atz>6 

We ascertain plausible scenarios by comparing the model results 
discussed in Sect. [3] with the derived values from observations 
for specific quantities of individual QSOs listed in Table [2] The 
calculated values for diverse properties such as M<j, M», Mh 2 , 
metallicity, and SFR from the models discussed below, which 
best match the QSOs, are listed in Table [3] The corresponding 
model parameters, and all models which match the discussed 
properties within the range defined by observations, are summa- 
rized in Table [4] 

We find that at an epoch of 30 Myr the models with an initial 
mass of the galaxy of M m \ = 1 x 10 11 M , an initial SFR of ij/ m \ 
= 3 x 10 3 M yr~' and either a Larson 2 IMF, a top-heavy or 
a mass-heavy IMF reproduce the observed quantities of some 
QSOs at z > 6 in the case of a 'maximum' SN efficiency. 

In particular, the model with a top-heavy IMF is best appli- 
cable to the QSO Jl 148+5251. The amount of dust reached is 
between 3.1-5.1 x 10 8 M for dust destruction in the ISM with 
Af c i = 100-0 M . A stellar mass of M* ~ 3.5 x 10 10 M is ob- 
tained. The metallicity in the system is ~ 2 Z and a SFR of 
~ 1600 M yr~' could be sustained. This model is also favored 
given its values of a an d rj g ^ 2 . The hig her H 2 mass of Mh, = 3.7 
x 10 10 M derived bv lRiechers et al.l d2009t) leads to rj gfi2 < 2 
and a ~ 1.4 M (K km s _1 pc 2 ) -1 . However, such a galaxy with 
Mi n i = 1 x 10 11 M implies that the dynamical mass is larger 
than the derived M dvn of ~ 5.5 x 10 10 M (for a i = 65°) by 
IWalter et all ((2004). While none of the model s for M iai = 5 x 
10 10 M , which was used bv lDwek et al.l (120071) . can be applied, 
a lower inclination angle similar to what has been adopted for 
the other QSOs might be considered. 

Another possible match with the properties of J 1148+5251 
is achieved by the same set of values for M; n i, i/^i, SN efficiency 
and IMF at an epoch of 100 Myr. The calculated stellar mass 
is within the estimated range from observations and the dust 
mass is ~ 2.4-8.9 x 10 8 M , depending on M i. However, the 
SFR dropped to ~ 1000 M yr _1 , while the metallicity increased 
to ~ 5 Z . In view of the lower SFR reached by these models 
than suggested by observations at epochs either 30 or 100 Myr, 
a higher initial SFR than the 3 x 10 3 M yr~' might be con- 
ceivable. In Fig. [3] one notices that a longer evolution with the 
same (or lower) initial SFR as used here does not lead to a better 
agreement with observations, since this results in an even lower 
SFR and higher metallicity. 

In view of this we find that this scenario at an epo ch of 100 
Myr is more appropriate for the QSOs J 1048+463 7 dFan et al.l 
2003) at z = 6.23 and J2054-0005 djiang et al.ll2008l) at z = 6.06. 
For the latter QSO a fine tuning of the epoch to 70 Myr results in 
a better match. At this epoch we obtain a SFR of 1 150 M yr~' 
and a metallicity of ~ 4.4 Z . The amount of dust is Mj ~ 2.7 x 
10 8 M (for M c i = 100 M ), while the stellar mass is M, ~ 4.7 
x 10 10 M . The lower derived L' CQ{l Q) leads to ^h, ~ 3—4 in 
case a = 0.8-1 M (K km s~' pc 2 ) -1 is applied, while for rj g h, ~ 
2 a value for a of ~ 1.6 would be required. For J1048+4637 the 
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Table 2. Observed properties of quasars at z > 6.0 



Object 


z 


^cor 


1-0) 


SFR 


M d 




M dyn sin 2 ( 


Ref. 






lO^Kkms-' pc 2 


M Q yr 1 


10 s 


M G 


10'° M Q 


10 10 M G 




Jl 148+5251 


6.42 


3.0 + 


0.3 


2380 


5.9: 


t 0.7 


2.4/3.7" 


4.5 


1,2,3,4 


J 1048+4637 


6.23 


1.2 ± 


0.2 


650 


4.3 : 


t0.6 


1.0 


4.5 


1,2,3 


J2054-0005 


6.06 


1.5 ± 


0.3 


1180 


3.4: 


h0.8 


1.2 


4.2 


5,2,3 


J0840+5624 


5.85 


3.2 + 


0.4 


1460 


4.7+ 0.9 


2.5 


24.2 


6,2,3 



References. (1 ) iFan et alj d2003h : (2) IWang et alj d2010h : (3) iMichalowski et alj |2010); (4) IWalter et alJ d2004h : (5) Uiang et alJ <2008h : (6) 
iFan et al.l J2006I) 

Notes. (a) Mh 2 = 3.7 x 10 10 M G deduced from [C I] line detections bv lRiechers et alj d2009h 



Table 3. Calculated properties from the best matching models of z > 6 QSOs from our sample 



Object" 


SFR 


M d 


M, 


z 




'/g.H, 


M H , 6 




M Q yr 1 


10 8 M Q 


10 10 M Q 


z Q 






10 10 M Q 


J 1148+5 25 1(A) 


1600 


3.1-5.1 


3.5 


2 


0.8-2.3 


3.0-1.0 


2.16-6.2 


J1148+5251(B) 


1000 


2.4-8.9 


5.4 


5 


0.8-1.55 


2.0-1.0 


2.10-4.1 


J1048+4637(A) 


1000 


2.4-8.9 


5.4 


5 


0.8-2.8 


3.4-1.0 


1.2—4.2 


J1048+4637(B) 


610 


3.5 


2.8 


3.4 


0.8-4.5 


5.8-1.0 


1.2-6.7 


J2054-0005 


1150 


2.7 


4.7 


4.4 


0.8-3.2 


3.0-1.0 


1.2-4.8 


J0840+5624(A) 


1500 


2.1 


11.0 


4 


0.8-4.6 


7.0-1.2 


2.5-14.7 


J0840+5624(B) 


1400 


4.8 


20.0 


5 


0.8-4.6 


10-1.8 


2.5-14.7 



Notes. (fl) All models are calculated for a top-heavy IMF. Capital letters in brackets (A,B) signify the different models (see corresponding models 
in Tabled for the same object. ( ' The ranges of rj g ,h 2 and Mh 2 corresponds to the range of a, which is between the commonly used value of a = 
0.8 and the possible upper limit. 



model for a lower initial SFR of 0^ = 10 3 M yr~' might be an 
option. The SFR is ~ 610 M Q yr~' and the metallicity is ~ 3.4 
Zq. While the stellar mass remains low, M, ~ 2.8 x 10 10 M Q , a 
dust mass of ~ 3.5 x 10 8 M is obtained for a 'maximum' SN 
efficiency and moderate dust destruction in the ISM. However, 
for a = 0.8-1 M (K km s _1 pc 2 ) 1 the gas-to-H 2 ratio is ~ 5-6, 
since for the lower initial SFR the system at this epoch is less 
exhausted. 

At either the same or a later epoch the more massive galaxies 
with M ini = 3-5 x 10 11 M , an initial SFR of iff m = 3 x 10 3 M 
yr _I and IMFs biased towards higher stellar masses are applica- 
ble to some z ~ 6 QSOs. The stellar mass, metallicity, and SFR 
of these systems are in agreement with observations, with either 
top-heavy IMFs or a mass-heavy IMF leading to the best results. 
The amount of dust can be produced by SNe with a 'high' SN 
efficiency and M c \ < 100 M , although the dust masses reached 
are at the lower limit. 

At an epoch of 170 Myr the system with M ini = 3 X 10 11 
M Q is plausible for the QSO J0840+5624 (IFan et alJl2006h at z 
= 5.85, if an inclination angle higher than the assumed 40° is 
assumed. The SFR is ~ 1500 M yr _I and the metallicity is ~ 
4 Z . The stellar mass is around 1.1 X 10 11 M . The amount of 
dust obtained with a 'high' SN efficiency is 2.1 X 10 8 M , while 
with the 'maximum' SN efficiency the dust mass exceeds a few 
times 10 9 M (as already at an epoch of 100 Myr). However, 
for a ^co(i-o) = 32 x 1qW k km S_I 

pc 2 as derived for this 
QSO the gas-to-H 2 ratio of i] gM , ~ 5-7 for a = 0.8-1 M (K 
km s -1 pc 2 ) -1 is higher than for the less massive galaxies. In 
case of a lower ^ g H 2 of ~ 2, a ~ 2.7 M (K km s pc 2 ) -1 is 
required. The larger galaxy with M m \ = 5 x 10 11 M , i/rj n ; = 3 x 
10 3 M yr -1 and top heavy IMF can account for the observed 
quantities at an epoch of 400 Myr. The amount of dust reached 



with a 'high' SN efficiency is ~ 4.8 x 10 8 M and the SFR is ~ 
1400 M yr -1 . The metallicity and stellar mass are in agreement, 
but the fraction of M H2 is around 1/10 for a = 0.8 M (K km s _1 
pc 2 )" 1 , while a ~ 4 M (K km s _1 pc 2 ) -1 is needed for of 
~ 2. A higher amount of Mh, as denoted by the higher value of 
a in these massive galaxies might be possible. For example, the 
presence of large amoun ts of cold and low-excited molecular gas 
have been suggested by Papadopo ulos et af] d2001l) for the QSO 
APM 08279+5255 at z = 3 .9 1 . 



4.2. SN efficiency and mass of the galaxy 

Our calculations show that with increasing M m \ (and fixed ifr m i, 
IMF) the SN dust production efficiencies can either be lowered 
or the degree of dust destruction increased in order to reach the 
required large dust masses. This is best demonstrated by models 
for the most massive galaxies with M m { = 3-13 x 10 11 M in 
which a 'high' SN efficiency is sufficient in case of moderate to 
no dust destruction. 

However, the largest system with M m \ = 1.3 x 10 12 M ex- 
ceeds the plausible dynamical masses derived from observations 
of QSOs at z > (5) 6 by more than an order of magnitude. 
Moreover, our computed models show that at least one of the 
properties of either SFR, Z or M« are not in agreement with ob- 
servations at any epoch fo r any as sumption of either the initial 
SFR or the IMF (see also Paper I). Additionally the values for 
77 g H 2 remain very high even for a = 4.6 M (K km s~' pc 2 ) -1 . 
W e therefore conclude t hat such a massive system as advocated 
bv lValiante et al.l (|2009), cannot be applied to QSOs at z > (5) 
6. Although systems with M; n ; = 3-5 x 10 11 M are appropriate 
for some QSOs at z < 6, such massive systems can only be ap- 
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Table 4. Models" which match the observed range of properties of z > 5 QSOs 



Epoch 


Mm 




SN efficiency 


M a 


IMF 


z > 6 QSOs from our sample 6 


M G 


10 3 M Q yr 1 




M 










30 Myr 


5 x 10 IU 


3 


max 


100 


top-heavy, Larson 1, 2, mass-heavy 








3 


max 





top-heavy, Larson 1, 2 








3 


max 





mass-heavy 








1 


max 


0,100 


top-heavy, Larson 2 






1 x 10" 


3 


max 


100 


Larson 2, mass-heavy 








3 


max 





mass-heavy 








3 


max 


100-0 


top-heavy 


Jl 148+525 1(A) 






3 


high 





Larson 2 




70 Myr 


1 x 10" 


1 


max 


100-0 


top-heavy 


J2054-0005 


100 Myr 


5 x 10'° 


3 


max 


100 


Larson 2 








3 


max 





top-heavy, Larson 2 








1 


max 


100 


Larson 2 






1 x 10" 


3 


max 


100 


Larson 2 








3 


high 





Larson 2 








3 


max 


100-0 


top-heavy 


Jl 148+525 1(B), J1048+4637(A) 






1 


max 


100 


Larson 2 








1 


max 


100-0 


top-heavy 


J1048+4637(B) 




3 x 10" 


3 


max 


800 


Larson 1 , mass-heavy 








3 


max 


800 


top-heavy, Larson 2 








3 


max 


100 


Larson 1 , mass-heavy 








3 


high 





mass-heavy 








3 


high 





top-heavy, Larson 2 








1 


max 


800 


Larson 2 








1 


high 





Larson 2 






5 x 10" 


3 


max 


800 


mass-heavy 








3 


max 


800 


top-heavy, Larson 2 








3 


high 


100 


top-heavy, mass-heavy 








3 


high 


100 


Larson 2 








3 


high 





mass-heavy 








3 


high 





top-heavy, Larson 2 




170 Myr 


3 x 10" 


1 


high 


100-0 


top-heavy 


J0840+5624(A) 


400 Myr 


5 x 10" 


1 


high 


100-0 


top-heavy 


J0840+5624(B) 



Notes. (n) All models which match the observed range of all properties of z > 6 QSOs are marked in bold face. <b) Capital letters in brackets signify 
that different models (A,B) are applicable for the same object 



plied to QSOs > 6 when the inclination angle is lower than the 
assumed average angle. 

The models which best reproduce the observed properties of 
QSOs > 6 are for a galaxy with M ini = 1 x 10 11 M G , but neces- 
sitate a 'maximum' SN efficiency and/or a moderate amount of 
dust destruction. The overall rapid evolution of dust and some 
properties in these models indicates that such QSOs could pos- 
sibly be present at a higher redshift than z > 6.4. An interest- 
ing example at a lower redshift of z = 1.135 is the ULIRG SST 
J 1604+4304, which shows proper ties similar to the considered 
high-z QSOs. iKawara et alj (1201 Oh reported a dust mass in this 
ULIRG of 1-2 x 10 8 M Q , a metallicity of around 2.5 Z and 
estimated the age of the stellar population to be 40-200 Myr. 

The possibility o f moder ate dust destruction in the ISM was 
already discussed in Paper I. We found that the amount of dust 
for most models better coincide with observations for M c \ < 
100 M , which would be in agreement with the values of M c \ 
of 50-70 M der ived for a multiphase ISM (e.g.. lMcKedll989t 
iDwek et al.ll2007h . 

The 'maximum' SN efficiency might be problematic. There 
is onl y little observational ev i dence that SN can be very efficient 
(e.g. IWilson & Batrlall2005t iDouvion et al.ll200lt iDunne et al.1 
2009), and theoretical models pred ict significant dust destruc- 
tion in reverse shocks of SNe (e.g.. Bianc hi & Schne ider 2007; 



iNozawa et all 120071 l2010h . On the other hand, these models 
also show that the effectiveness of dust destruction depends 
on various properties such as the geometry of the shocks, the 
density of the ejecta and the ISM, the size and shape of the 
grains, clumping in the SNe ejecta, and different SN types. 
In addition there is some observational evidence that Type Iln 
SNe and sources such as l uminous blue variables are possi- 
bly efficient dust p roducers dFox et al. 2009} ISmith et al.l l2009: 
Gomez et al. 2010). While dust production and destruction in 
SNe is yet unresolved, a 'maximum' SN efficiency cannot be 
ruled out (e.g., Gall et al. in prep). Alternatively, either dust 
formation in the outflowing wind s of QSOs or gra i n growth in 
the IS M migh t be a n option (e.g., Elvis et al.|l2002t IDwek et al.l 
| 2007tlDrainell2009l:lMichaIowski et al.ll2010t IPipino et alj|201 \i 
IDwek & Cherchnefl 1201 Oh as supplementary or primary dust 
sources. However it remains to be investigated, if dust grain 
growth can be as efficient as required under the prevailing con- 
ditions of high star formation activity and a short time span. 
Typical grain growth timescales in molecular clouds are of or- 
der 10 7 yr, but depending on the density and metallicity these 
can possibly be sh orter (e.g., Hirashita 2000; Zhukovsk a et al.l 
2008; Draine 2009). The fact that the starburst is assumed to oc- 
cur in an initially dust free galaxy implies that heavy elements 
first need to be ejected into the ISM before grain growth can take 
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place. In forthcoming work we will further develop the model to 
investigate the impact of different infall and outflow scenarios 
on the evolution of the amount of dust and various properties of 
a galaxy. 
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